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We construct a nonequilibrium theory for the charge transfer through a diffusive array of alter- 
nating normal (N) and superconducting (S) islands comprising an SNSNS junction, with the size 
of the central S-island being smaller than the energy relaxation length. We demonstrate that in 
the nonequilibrium regime the central island acts as Andreev retransmitter with the Andreev con- 
versions at both NS interfaces of the central island correlated via over-the-gap transmission and 
Andreev reflection. This results in a synchronized transmission at certain resonant voltages which 
can be experimentally observed as a sequence of spikes in the differential conductivity. 
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An array of alternating superconductor (S) - normal 
metal (N) islands is a fundamental laboratory represent- 
ing a wealth of physical systems ranging from Joseph- 
son junction networks and layered high temperature su- 
perconductors to disordered superconducting films in the 
vicinity of the superconductor-insulator transition. Elec- 
tronic transport in these systems is mediated by An- 
dreev conversion of a supercurrent into a current of 
quasiparticles and vice versa at interfaces between the 
superconducting and normal regions A fascinat- 

ing phenomenon benchmarking this mechanism is the 
enhancement of the conductivity observed in a single 
SNS junction at matching voltages constituting an in- 
teger (m) fraction of the superconducting gap, V = 
2 A/ (em) ^ due to the effect of multiple Andreev reflec- 
tion (MAR) [||. The current- voltage characteristics of 
diffusive SNS junctions were studied in detail in [|[ ||. 
Recent experimental findings [||-D posed the question 
about the transport properties of large arrays comprised 
of many SNS junctions. In this Letter we develop a 
nonequilibrium theory for the current- voltage character- 
istics of a series of two diffusive SNS junctions, i.e. an 
SNSNS junction. The normal parts have, in general, dif- 
ferent resistances and are coupled via a small supercon- 
ducting granule, S^ . The focal point is the construction 
of a nonequilibrium circuit theory for the charge transfer 
across the superconducting island. We demonstrate that 
the wedging of S^ into the normal part of an SNS junction 
leads to nontrivial physics and the appearance of a new 
distinct resonant mechanism for the current transfer: the 
Synchronized Andreev Transmission (SAT). In the SAT 
regime the processes of Andreev conversion at the bound- 
aries of the central superconducting island are correlated: 
as a quasiparticle with energy e hits one NS^ interface, 
a quasiparticle with the same energy emerges from the 
other SpN interface into the bulk of the normal island 
(and vice versa, see Fig. 1). This energy synchroniza- 
tion is achieved via over-the-gap Andreev processes 0, 



which couple the MAR processes occurring in each of the 
normal islands and make the quasiparticle distributions 
at the central island essentially in nonequilibrium. Effec- 
tiveness of the synchronization is controlled by the energy 
relaxation lengths of both, the quasiparticles crossing S,, 
with energies above A, and of quasiparticles experiencing 

N, 




FIG. 1: (color online). Diagrams of the SAT processes for 
n = 1 in Eq. (|l|) in an SNSNS junction with normal resis- 
tance ratio Ri/R'z =3/4 (depicted through the 3/4 ratio of 
the respective lengths of the normal regions). The thick solid 
lines represent the synchronized quasiparticle paths connect- 
ing the points of singularity in the density of states corre- 
sponding to energies e = ±A at the electrodes S^ and Sj, . 
The paths ABCD (upper panel) and D'C'B'A' (lower panel) 
correspond to the electron- and hole trajectories, respectively. 
Synchronization of the energies of the incident and emitted 
quasiparticles at points B and C (B' and C') is shown by 
arrows. SAT is realized by trajectories passing through the 
singular points ±A of the central island S,-, and including 
above-the-gap Andreev reflections. Trajectories synchroniz- 
ing other transmissions across S(-. and those of higher orders 
are not shown. Note, that voltage drops eVi = 6A/7 and 
eV2 = 8A/7 are not MAR matching voltages of either the 
S^NiSc or the ScN2Spt part. 
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MAR in the normal parts. The SAT processes result in 
spikes in the differential conductivity of the SNSNS cir- 
cuit, which appear at resonant values of the total applied 
voltage Vtot defined by the condition 
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with integer n, irrespectively of the details of the distribu- 
tion of the partial voltages at the two normal islands. As 
we show below, the SAT-induced features become domi- 
nant in large arrays consisting of many SNS junctions. 

The model. — We consider the charge transfer across an 
SlNiS^N2Sj^ junction, where S^, S^, and Sj, are super- 
conducting islands with identical gap A. We assume the 
size of the central island to be much larger than the 
superconducting coherence length ^, hence processes of 
subgap elastic cotunneling and/or direct Andreev tunnel- 
ing Q are not effective. In general this condition ensures 
that is large enough so that charges do not accumu- 
late in the central island and Coulomb blockade effects 
are irrelevant for the quasiparticle transport. At the same 
time is assumed to be less than the charge imbalance 
length, such that we can neglect the coordinate depen- 
dence of the quasiparticle distribution functions across 
the island S^ . Additionally, the condition £^ ^ L^, 
where £1; is the energy relaxation length, implies that 
quasiparticles with energies e > A traverse the central 
superconducting island S^, without a loss of energy. The 
normal parts Ni and N2 are the diffusive normal met- 
als of length Li^2 > ^, and Li_2 > L,^, L,^ = ^JhD^Je, 
where is the diffusion coefficient in the normal metal. 
We assume the Thouless energy, iJxh = ^^n/^i.2: to 
be small, i^Th ^ A, and not to exceed the character- 
istic voltage drops, E'xh < eVi.2. These conditions are 
referred to as incoherent regime where, in particular, 
the Josephson coupling between the superconducting is- 
lands is suppressed. We let the energy relaxation length 
in the normal parts Ni and N2 be much larger than their 
sizes, thus quasiparticles may experience many incoher- 
ent Andreev reflections inside the normal regions. 

The current transfer across the SNSNS junction is de- 
scribed by quasiclassical Larkin-Ovchinnikov (LO) equa- 
tions for the dirty limit |10[: 



-i[i7cffO,G] = VJ, J-n: 
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where iJoff' = ^{i<yzdt — fcro + A), J = DG o VG is 
the matrix current, the subscripts "S" and "N" stand 
for superconducting and normal materials, respectively, 
"o" is the time-convolution, ct" (a = {0,1,2,3}) are the 
Pauli matrices, A = ia-^lmA + io-y ReA, and R is the 
resistance of an NS interface. The diffusion coefficient D 
assumes the value Z)„ in the normal metal and the value 
in the superconductor, and ip is the electrical poten- 
tial which we calculate self-consistently. The unit vector 



n is normal to the NS interface and is assumed to be 
directed from N to S. The momentum averaged Green's 
functions G(r, t, t') are 2x2 matrices in a Keldysh space. 
Each element of the Keldysh matrix, labelled with a hat 
sign, is, in its turn, a 2 x 2 matrix in the electron-hole 
space: 



G = 
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r is the spatial position and t, t' arc the two time ar- 
guments. The Keldysh component of the Green's func- 
tion is parametrized as G^ = G^ o / — / o G^, 
where / is the distribution function matrix, diagonal in 
Nambu space, / = diag [1 — 2nc, 1 — 2nh], n^ih) is the elec- 
tron (hole) distribution function. In equilibrium rio(h) 
becomes the Fermi function. And, finally, the Green's 
function satisfies the normalization condition & = 1. 

The edge conditions closing Eqs. (^) are given by the 
expressions for the Green's functions in the bulk of the 
left (L) and right (R) superconducting leads: 



G 



L(R) 



(t, t') = e-*^L(R)tf3/fiQ^(^ _ i')g^/XL(R)t'f3/S ^ 



the chemical potentials are /i^ = and /ip^ = eV . Here, 
Go(i) is the equilibrium bulk BCS Green's function. 

The current density is expressed through the Keldysh 
component of J as 

I(t,r) = ^TrajK(t,t;r) j de[I,{e) + h{e)] , 

(4) 

where the spectral currents and /h are the time 
Wigner-transforms of top- and bottom diagonal elements 
of the matrix current 3^^\ representing electron and hole 
quasiparticle currents, respectively. In the bulk of a nor- 
mal metal /e = a^Vrie and /h = cr^Vnh. 

The distribution functions of quasiparticles at the cen- 
tral island S,, are essentially is nonequilibrium. To take 
this into account we define quasiparticle spectral cur- 
rents at NS interfaces by the Keldysh component of the 
boundary conditions for Eqs. (0). These nonequilibrium 
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FIG. 2: (color online). Effective circuit representing current 
conversion at the interfaces of the central superconducting 
island S^. Resistors, i?p and i?^ stand for an Andreev- and 
a normal processes, respectively. 
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boundary conditions have a form of Kirchhoff 's laws for 
the circuit shown in Fig.|2[ The electron and hole distri- 
bution functions take the role of voltages at the nodes. 
For illustration we write down the equation for an elec- 
tronic spectral current flowing into the lower left corner 
node (Kirchhoff's laws at the other corner nodes have a 
similar form): 



1 



+ 



Rr\e + u) 



4''(£ + 2M)-nW(e) . (5) 



The interjacent resistances, -Rq(p), are defined as 
R~l^^{£) = ± ^+^(e)}/2 where R±{e) are the 

special functions characterizing transparencies of the in- 
terfaces and tabulated in Q]. At high energies, \e\ ^ A, 
R+{e) R-is)', at small energies, |£| ^ A, R-{e) di- 
verges, while i?+ (e) — i?_ (e) remains finite; and there are 
singularities in deR±{e) at |e| = A, which are of the same 
origin as those in the BCS density of states and reflect 
the fact that quasiparticles cannot penetrate the super- 
conductor below the gap. The "bars" indicate that the 
respective resistances and the distribution functions are 
renormalized by the proximity effect. 

To derive the current-voltage characteristics for the 
general case of an asymmetric nonequilibrium SNSNS 
junction with different resistances of the normal regions, 
we construct a nonequilibrium circuit theory allowing for 
an analytical solution of the nonlinear nonuniform matrix 
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FIG. 3: (color online). Effective circuit for an SNSNS junc- 
tion. To simplify the graph, resistors R^^^-f ~ (l/2){[^+"^^ + 

-R+'^'l ""^ ± [^R^i"^' + .R*^"^']^^} representing Andreev and nor- 
mal processes at the NS interfaces of the central island are 
introduced. 



Eqs. (g) for Keldysh-Nambu Green's functions. The di- 
agrammatic mapping of Eqs. is realized by an equiv- 
alent circuit shown in Fig.^. The Kirchhoff's equations 
for the potential distribution in the circuit of Fig. ^, give 
the recurrent relations: 



TZie, -u, -V)hie) - p^°\e - n)h{e - 2u) 



{e)I,{e)~ p(^\e~V)h{e-V) (6) 



7^(e, u, F)/e (e) - (e + u)/h (e + 2u) 
-p(^)(e)/h(e) - p^^\e + V)Ue + 2V) 



(7) 



where the electric potential u of the Sc island is cal- 
culated self-consistently from the electroneutrality con- 
dition, u — (7r/8)TrG"^. The effective resistance is 



n = Ri+R2 
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where p 
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(l/2){i?^'''^ - iiL'"'^}. Solutions of Eqs.(|) yield the re- 
quired I-V characteristics for an asymmetric SNSNS 
junction. To verify our formulas we note that at large 
quasiparticle energies, |e| » A, the total resistance "R. 
reduces to the normal resistance of the array, whereas 
p^"') and p^'^') vanish. Then we find from Eqs.(||)-(0) that 
/h(e) - {n^(e)-n^(e-V)\l'R, and /e(£) = {n^(e + V)- 
np(e)]/7?., which together with Eq.(^) reproduce Ohm's 
law, I = VjH,. The constructed diagram, Fig. ^, is an el- 
emental building unit for a general nonequilibrium quan- 
titative theory of SNS arrays comprised of many SNS 
junctions. 

Results and discussion. — The calculation of the 
current-voltage characteristics I(V) requires the numer- 
ical solution of the recurrent relations, Eqs. (^)-(0). To 
this end, we have developed a computational scheme al- 
lowing to bypass instabilities caused by the non-analytic 
behavior of the spectral currents /e(h)(£). We first fix 
some chosen energy e, identify the set of energies con- 
nected through the equations in the given energy inter- 
val, and solve the resulting subsystem of equations. We 
then repeat the procedure, until the required energy res- 
olution of Se = 10~^A is achieved. Typically, up to 10^ 
linear equations had to be solved for every given voltage, 
but the complexity of the coupled subsystem depends on 
the commensurability of u and V . 

Figure ^ shows the comparative results for the SNSNS 
junction and two SNS junctions in series. The latter cor- 
responds to the case where the size of the central island 
well exceeds the energy relaxation length, Lc > le- We 
display the differential resistances as functions of the ap- 
plied voltage, which demonstrate the singularities in An- 
dreev transmission more profoundly than the I-V curves. 
There is a pronounced SAT spike in the dV/dl for an 
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FIG. 4: (color online). Left panel: Differential resistances as 
functions of the applied voltage Vtot (around n = 1 in Eq. (1)) 
for the SN1SN2S junction. The fractions 3/4, and 4/5, and 
9/10 represent the ratios of resistances of the normal regions, 
Ri/R2- dV/dl of SN1SN2S junction demonstrates the pro- 
nounced SAT spike at Vtot ~ 2A/e, irrespectively of the par- 
tial voltage drops. The SAT spike is sandwiched between the 
two additional spikes corresponding to individual MAR pro- 
cesses occurring at junctions SNiS and SN2S for mi,m2 = 2. 
The voltage positions of these features depends on R1/R2. 
Right panel: The corresponding dV/dI{Vi + V2) for two SNiS 
and SN2S junctions in series as they would have appeared in 
absence of the synchronization process, i.e. when Lc > £e- 
These dV/dl were calculated following (with transmissiv- 
ity W=l). 



SNSNS junction at Vtot — 2A/e. The spike appears ir- 
respectively of the partial voltage drops in the normal 
regions and is absent in the corresponding curves repre- 
senting two individual MAR processes at the junctions 
SNiS and SN2S. 

The resonant voltages of the SAT singularities can be 
found from the consideration of the quasiparticle trajec- 
tories in the space-energy diagrams. Such a diagram for 
the first subharmonic, n = 1 and ratio R1/R2 ~ 3/4 
is given in Fig.|^. A quasiparticle starts from the left 
superconducting electrode with the energy e = —A to 
traverse Ni, and the quasiparticle that starts from the 
central island Sc with the same energy as the incident 
one to take up upon the current across the island N2, 
and hit Sr with the energy e = A (the ABCD path, the 
corresponding path for the hole is D'C'B'A'). In gen- 
eral, relevant trajectories yielding resonant voltages of 
Eq. (|^) have the following structure: they start and end 
at the BCS quasiparticle density of states singular points 
(e ~ ±A), contain the closed polygonal path, which in- 
clude MAR staircases in the normal parts and over-the- 



gap transmissions and Andreev reflections, and pass the 
density of states singular points at the central island. 
Apart from the main singularities [Eq. (|l|)], additional 
SAT satellite spikes appear at V = {2A/e){p + q)/n, 
where p/q is the irreducible rational approximation of 
the real number r = R1/R2, (we take i?i < R2), and 

" > {p + q)- 

The achieved qualitative understanding enables us to 
observe that the manifestations of the SAT mechanism 
in an experimental situation becomes even more pro- 
nounced with the growth of the number of SNS junctions 
in the system. To see this, let us assume that the resis- 
tances of the normal islands in a chain of SNS junctions 
are randomly scattered around their average value Ro 
and follow Gaussian statistics with the standard devia- 
tion ctr = cri?0: where a is dimensionless. Accordingly, 
the dispersion of the distribution of the MAR resonant 
voltages is characterized by the same a, and the MAR 
features get smeared. Let us distribute the voltage drop 
2A/e among the n successive islands. Then the quasi- 
particle SAT path starts at the lower edge of the super- 
conducting gap at island j, traverses n — 1 intermediate 
superconducting islands and hits the edge of the gap at 
the j + n-th island in the chain. The standard deviation 
of the voltage drop on the n islands grows as ^/n result- 
ing in a voltage deviation per one island oc 1/y/n, i.e. the 
dispersion of the distribution of Vn drops with increasing 
n: (TsAT = o-/y/n- In contrast to the MAR-induced fea- 
tures, with an increase of n, the subharmonic spikes at 
voltages Vn per junction due to SAT processes become 
more sharp and pronounced. 

In conclusion, we have developed a nonequilibrium the- 
ory of charge transfer across a superconducting island 
and found that the island acts as Andreev retransmit- 
ter. We have shown that the nonequilibrium transport 
through an SNSNS array is governed by synchronized An- 
dreev transmission with correlated conversion processes 
at the NS interfaces. The constructed theory is a funda- 
mental building unit for a general quantitative descrip- 
tion of a large array consisting of many SNS junctions. 
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